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Abstract. We consider the long-time behaviour of a branching random walk in ran- 
dom environment on the lattice The migration of particles proceeds according to 
simple random walk in continuous time, while the medium is given as a random poten- 
tial of spatially dependent killing/branching rates. The main objects of our interest 
are the annealed moments {mn), i.e., the p-th moments over the medium of the n-th 
moment over the migration and killing/branching, of the local and global population 
sizes. For n = 1, this is well-understood jGM98j . as mi is closely connected with the 
parabolic Anderson model. For some special distributions, jABMYOO] extended this 
to n > 2, but only as to the first term of the asymptotics, using (a recursive version 
of) a Feynman-Kac formula for m n . 

In this work we derive also the second term of the asymptotics, for a much larger 
class of distributions. In particular, we show that (mS) and {m r [ p ) are asymptotically 
equal, up to an error e°(*'. The cornerstone of our method is a direct Feynman- 
Kac- type formula for m n , which we establish using the spine techniques developed in 
|HR12| . 

1. Introduction 

Random processes in random surroundings are under investigation for decades. Examples of such 
processes include (1) trajectories of random walks and Brownian motion in random environment 
with the focus on laws of large numbers and central limit theorems or even invariance principles, (2) 
heat equation and other partial differential equation systems in random potential with the focus on 
intermittent behaviour, (3) polymer measures and directed percolation in random medium with the 
focus on free energies. In all these models, a rich phenomenology of asymptotic behaviours arises, 
which is not shared by the original model in non-random (homogeneous) surrounding, and in most of 
these models the research goes on, as many of the main features have not yet been properly understood. 

Another large class of random processes that develop striking properties in random surroundings 
is the class of branching processes. Let us briefly describe some of the work that was carried out about 
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these models. Branching discrete random walks on Z, d with time-space i.i.d. offspring distributions 
were studied in the context of survival properties, global/local growth rates and diffusivity; and their 
connections to the directed polymers in random environment, see e.g. |BGK05( \Y08\ ICYllj . Detailed 
analyses of recurrence/transience properties of discrete-time branching Markov chains with only space- 
dependent environment, which does not exhibit in general the the usual dichotomy valid for irreducible 
Markov chains, were carried out in [CMP981 iMPOOl IMP031 ICP07l lM08l IBGK091 [GMPVlO] . to mention 
some. The main techniques in these studies relate these models to the better-known random walk in 
random environments, using the spectral properties of underlying Markov process and studying the 
embedded Galton- Watson processes in random environment. 

In this paper, we study a branching random walk in random environment (BRWRE), where 
the particles move around in space like independent random walks in continuous time, and the 
killing/branching takes place in sites with a random site-dependent rate. We are interested in the 
long-time asymptotics of the annealed moments of any order of the local and global population sizes. 
As was explained in |GM 90] for the case of first moments, this question stands in a close connection 
with the description of the intermittent behaviour of the main particle flow, i.e., its concentration be- 
haviour in small islands. According to the best of our knowledge, this question for the higher moments 
has hardly been investigated for this model yet, the only example being jABMYOO] . In that paper, a 
deep relation between the moments of the BRWRE and the parabolic Anderson model is revealed and 
employed in order to analyse the annealed moments of the BRWRE, i.e., the p-th moments over the 
medium of the n-th moment over the killing/branching and migration of the total and local population 
size. It is the aim of the present paper to significantly increase the validity and the deepness of the 
results of [ABMYOOJ and to reveal the general mechanism that leads to the moment asymptotics. In 
contrast with |ABMY00| . we will be using probabilistic methods rather than PDE methods. 

1.1 Branching random walk in random environment 

Let us describe the model in more detail. The branching random environment on the lattice Z rf is 
a pair E = (^0^2) of two independent i.i.d. fields £0 = (^,o{y))y<=z d an d £2 = {S,2(y)) y e'E d °f positive 
numbers. Indeed, £o(y) and £2(2/) play the role of the rate of the replacement of a particle at y € Z d 
with or 2 particles, respectively. For n = 0, this is a killing, for n = 2, this is a binary splitting. (See 
Section [1.41 for more general branching mechanisms.) 

The probability measure corresponding to S is denoted Prob; expectation with respect to Prob 
will be written with angular brackets (•). For a given realization of S, the branching process with rate 
field S is now defined by determining that any particle located at a lattice site y € Z d , is subject to 
the killing/branching defined by the rates £0(2/) an d £2(2/)) and additionally each particle performs a 
continuous-time random walk on 7L d with generator re A, where re > is a parameter, and 

A/(x) = ]T [f(y) - f(x)] , for x € Z d , f G £ 2 (Z d ), 

is the standard lattice Laplacian. We write expectation with respect to a random walk with generator 
reA starting from x as ¥ x with corresponding expectation E z . We consider a localised initial condition, 
i.e., at time t = 0, there is a single particle at some site x € Z d . Probability and expectation w.r.t. the 
migration, branching and killing of the BRWRE are denoted by P x and E x , respectively, for fixed 
medium S. 

The description of the dynamics of the population is as follows. If a particle is at some time at 
some site y, then during a small time interval of length h, with probability re/i + o(h) it moves to a 
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neighbouring site chosen uniformly at random, with probability £2(2/)^ + o(h) it dies and is replaced 
by two descendant particles, and with probability £o(y)h + o(h) it is killed without producing any 
offspring. Finally, with probability 1 — (« + £2(2/) + £0(2/))^ + o{h) the particle experiences no changes 
during the whole time interval of length h. 

Let r](t,y) be the number of particles at time t G [0, 00) at y G Z d , and let n(t) = Yl y ez d vi^^v) 
be the total population size at time t. The main objects of interest in this paper are the quenched 
moments 

m n (t,x,y) = E x [i](t,y) n ] and m n (t, x) = E x [ V (t) n ], n G N, (1.1) 

i.e., the expected n-th powers of the local and global particle numbers, where the expectation is taken 
only over the migration and the killing/branching, for frozen killing/branching rates S. Note that, 
for n = 1, mi(i, x) is equal to the sum of m\{t,x,y) over y G Z d , but such a relation is not valid for 
n > 2. 

It will be the main purpose of the present paper to analyse the large-t asymptotics of the p-th 
moments of m n (t,x) and of m n (t,x,y), taken over the medium E. 

1.2 Connection with the parabolic Anderson model 

It is a fundamental knowledge in the theory of branching processes that the expected particle number 
satisfies certain partial differential equation systems. In our case, the characteristic system reads as 
follows. Put 

£ = £ 2 -£o, 

and fix y G then, (under certain integrability conditions, see |GM90j ) for fixed localised initial 
condition mi(0, •, y) = S y (-), the map (t, x) \-t mi(t, x, y) is the unique positive solution to the Cauchy 
problem for the heat equation with potential £, i.e., 

d 

—mi(t,x,y) = «;Ami(i, x, y) + £(x)mi(t, x, y), for (t, x) G (0, 00) x Z d . (1.2) 

Similarly, the map (t,x) h-> mi(t,x) is the unique positive solution of (jl.2p with delocalized initial 
condition mi(t, •) = !■■ 

The interesting feature in our case is that the potential £ is random, and here (jl.2p is often called 
the parabolic Anderson model. In fact, the operator kA+£ appearing on the right-hand side is called the 
Anderson operator; its spectral properties are well-studied in mathematical physics. Equation (ll.2p 
describes a random mass transport through a random field of sinks and sources, corresponding to 
lattice points z with £(z) < and £(z) > 0, respectively. We refer the reader to |GM90| . |M94j 
and |CM94| for more background and to |GK05| for a survey on mathematical results. We see two 
competing effects: the diffusion mechanism (Laplacian) tends to make the field mi flat, and the local 
growth (potential) tries to make it irregular. 

Furthermore, it is also widely known since long [GM90] that mi admits a representation in terms 
of the Feynman-Kac formula: 

m 1 (t,x,y)=E x [exv{J^(X s ) ds}6 y (X t )\, (t, y) G [0, 00) x Z d , (1.3) 

and the same formula without the last indicator for mi(t,x), where (^ s ) s g[o,oo) denotes a simple 
random walk with generator kA. Note that mi depends only on the difference £ of £2 and £o- 
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The asymptotics of the moments of mi were analysed in |GM98j for the interesting special case 
that the distribution of £ lies in the vicinity of the so-called double- exponential distribution with 
parameter p € (0, oo), 

Prob(£(z) > r) = exp{-e r ^}, r € (0, oo). (1.4) 

The precise assumption on £ can be written down in terms of the logarithmic moment generating 
function 

H(t) = log(e t ^), (1.5) 
which is assumed to be finite for any t > 0. 

Assumption 1.1. There exists p £ [0, oo] such that 

lim H ( ct )~ cH ( t ) =pc lo gc , 6(0,1). (1.6) 

t— s-oo f 

Under this assumption, it is proven in |GM98j that, for any x £ Z rf , as t -> oo, 

( m {(t, x)) = e H{pt) e -^x( P /^)pt+o(t) ^ p € Nj (L7) 

where x is defined as 

X {p) = \ mi[S{p)+pl{p)\. (1.8) 

Here "P(Z) denotes the space of probability measures on Z, and the functionals S,I: "P(Z) — > M + are 
given by 

S(p) = Y,W^ x + l )-y^) 2 and X(/,) = -^Mx)log/.(x). (1.9) 

We have < xip) < 1 for p € (0, oo) and x(0) = and lmip^oo x(p) = x(°°) = 1- The right-hand 
side of (|1.7|) is also equal to the moments of m±(t,x,y) for any fixed x, y £ Z rf , as is seen from an 
inspection of the proof (see Remark 1.3 in [GM98] ). Also note that H{t) 3> 2dKx(p/K)pt for large t, 
that is, asymptotically the first term on the right-hand side of (jl.7p is much larger than the second 
term. 

Observe that the p-th moments of mi at time t behave like the first moment at time tp, up to 
the precision of (jl.7p . This can be easily guessed from a standard eigenvalue expansion for m\(t,x) 
in terms of the eigenvalues and eigenfunctions of kA in large t-dependent boxes with zero or periodic 
boundary condition; in fact, mi(t, x) is roughly equal to e tXl ^\ where X\(t) is the principal one. Then, 
obviously, m^(t,x) is roughly equal to e* pAl (*). 

1.3 Moments of the BRWRE 

Let us now turn to the main object of the present paper, the moments of m n for n > 2. We can 
formulate our main result. Recall our assumptions from the beginning of Section 11.11 We will also 
suppose that the branching rate £2(0) satisfies Assumption 11.11 In the case p = 00, we will need an 
extra assumption to avoid too large a growth of H2(t): 

Assumption 1.2. For any k € N, 

<6(0) fc e^) < <e^° V W as t -» 00. (1.10) 
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Theorem 1.3 (Moments of the BRWRE). Suppose that the logarithmic moment generating function 
H2 0/^2(0) satisfies Assumption li.il and, in the case p = 00, £2 also satisfies Assumption Fix 
x € Z d , the starting site of the branching process. Then, for any p, n € N, as t — > 00, 

{m p n {t, x)) = exp (H{npt) - 2dKx(p/n)npt + o(tj\ . (1.11) 
The same asymptotics holds true for (mn(t,x,y)) for any y € 

Note that the logarithmic moment generating function Hq of — £o(0) nas asymptotics -zHo(t) —> 
— essinf (£o(0)) £ (—00, 0] as t — > 00. Therefore, by independence of £2 and £oi the logarithmic moment 
generating function H of £2(0) — £o(0) also satisfies Assumption 11.11 an d this is crucial for the validity 
of Theorem 11.31 

In particular, Theorem 11.31 says that the p-th moments of m n at time t are equal to the first 
moment of mi at time tpn, up to the precision in (jl.llj) . i.e., 

(m p n (t,x)} = (m7(t,x))e° W = ( mi (tnp, x))e° W , t -»• 00. (1.12) 

This fact is not so easy to understand as for the case n = 1, see above. However, see Section [1.51 
for some heuristic remarks. The proof of Theorem 11.31 is in Section [3j 

The main tool of our proof is a Feynman-Kac-type formula for m n , which we will derive in 
Section see Theorem 12. 1L We are going to use probabilistic tools from the theory of branching 
processes, the main input coming from the many-to-few lemma of [HR12j . 

In [ABMYOO] there was a weaker version of (jl.lip derived; actually only the first term Q H i n P t ) ^ 
and this only for the rather restricted case of a Weibull distribution, where H{t) ~ Ct a for some 
C € (0,oo) and a G (l,co), a subcase contained in Assumption 11.11 in p = 00. On the other hand, 
they drop the assumption of independence and only assume spatial homogeneity of S. However, this 
result does not explain the spatial structure of the peaks of the moments of the population size, an 
information that is contained in the second term, as was discussed at length in |GM98j . The proof in 
[ABMYOOj is based on the fact that m n is the solution to an inhomogeneous Cauchy problem, where 
the inhomogeneity is a linear combination of products of mi, . . . , m n _i. Furthermore, they derived 
from this a Feynman-Kac formula for m n , which depends on that inhomogeneity and is therefore 
of recursive type. This made it rather difficult to identify the second term of the asymptotics. In 
contrast, we first derive a direct version of a Feynman-Kac-type formula in Theorem 12.11 and are then 
able to find the logarithmic asymptotics of the moments in much higher precision. 

1.4 More general branching 

Our Theorem 11.31 is formulated only for the special case of binary branching, but it can straightfor- 
wardly be extended to more general branching mechanisms, subject to additional conditions. Indeed, 
assume that the branching random environment is a family S = (£fc)fceN of i.i.d. fields = {£,k(y))y£z<i 
of positive numbers. Then £fc(y) is the rate for replacement of a particle at y by precisely k new par- 
ticles, i.e., a splitting into k particles. To exclude trivialities, we put £i(y) = f° r an Y V- The family 
(6c)fceN is n °t assumed to be i.i.d. Indeed, we at least have to assume that the field 

00 

e(y) = E( fc - 1 )^(y) ( L13 ) 

fc=0 

is well-defined (i.e., absolutely convergent) almost surely. One possible choice could be = £pk with 
some probability distribution (pfc)fceN an d some positive i.i.d. field £. 
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Then, under the assumption that J2keN k n £,k(y) < °° almost surely (e.g., if = for all suffi- 
ciently large k), our Feynman-Kac-type formula for m n in Theorem 12.1 1 below extends to this more gen- 
eral setting, see Remark 12 .21 Furthermore, under suitable conditions on the moments of XlfceN ^™£fc(2/)> 
also the proof of Theorem 11.31 in Section [3] can be easily extended to this situation. In order to avoid 
cumbersome formulas, we abstained from writing down the details. 

1.5 Discussion 

Let us explain why the moment asymptotics of mn(t, x) are equal to the ones of m^ n (t, x), see fjX. 12|) . 
We do this for n = 2 and p = 1. Note that, according to Theorem 12.11 below, = mi + 777,2, where 

m 2 (t,x) = j\ x [e W {J\(X r ) dr + J t(X' r ) dr + j\{K) dr}2&(X a )] da, (1.14) 

where (X r ) r& y Qs \ an d (^' r )re[s,t] an d (^rOrektl are independent simple random walks, given X s , with 
generator At A, starting at X$ = x, and X' s = X'J = X s . In other words, these three random walks 
constitute a branching random walk with precisely one splitting at time s. The first part in the 
decomposition 7772 = 7771 + 777,2, corresponds to absence of splitting, and the second one to precisely 
one splitting. 

Let us consider the behaviour of the moments of 777 2 as t — > 00. The first observation is that the 
term 2£ 2 (Xf) should have hardly any influence. This is due to Assumption 11.21 which rules out cases 
of extreme growth of H2 (t) . One can expect from (|1.7p that the leading term of the expectation on the 
right-hand side of (|1.14p should be Q H ( 2t ~ s ) f corresponding to the total time s + (t — s) + (t — s) that 
the three random walks spend in the random environment. Since H(t) —¥ 00, this is clearly maximal 
for s ~ 0. Hence, the Laplace method gives that the main contribution comes from s ~ 0. Hence, 
the contribution comes mainly from a product of expectations over two i.i.d. copies (X' r ) re i £ and 
(X") re \Q )t ], i.e., from a term mf(t,x). 

In other words, it is favourable for the branching random walk to split as soon as possible into 
two copies and to travel through the environment with these copies for a long time. The deeper reason 
for this is that the potential £ assumes extremely high values in some part of the space, where the two 
copies collect much of them. This effect seems to be present as soon as esssup(£(0)) is positive, and 
it should be turned into its opposite if esssup (£(0)) is negative. More precisely, for such potentials, 
we expect that (m n (t, x)) w (mi(t,x)). We expect that, for all four classes of potentials in the 
classification made in [HKM06] , a version of Theorem 11.31 can be deduced from Theorem 12.11 



2. FEYNMAN-KAC-TYPE FORMULA FOR 777 n VIA SPINE TECHNIQUES 

In this section, we derive a Feynman-Kac-type formula for 777 n , almost surely with respect to the 
branching rates £2 and killing rates £()■ Ovrr main result of this section appears in Theorem 12.11 below. 
We will use the spine techniques of [HR12j . This requires the introduction of a branching random 
walk (BRW) in 7L d with time interval [0, t] with up to 77 — 1 splitting events. In order to express this 
BRW, we will need the following ingredients. 

(i) a tree that expresses the branching structure, 

(ii) an ordering of the splitting sites of the tree to express their order in time, 

(iii) a time duration attached to each bond, 

(iv) an expectation over a simple random walk bridge attached to each bond. 
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In order to keep the notation simple, we restrict to binary branching. See Remark 12 .21 for more general 
branching mechanisms. 

We need some notation from the theory of trees. Let G = (V, E) be a finite graph with V the 
set of vertices and E the set of edges. G is a tree if it is simple, connected and has no cycles. Let 
us assume that G is a rooted tree, i.e., a tree with a root G V . This induces a natural ordering of 
vertices, namely, for u, v G V we say that u ■< v if the unique path from the root to v contains u. In 
particular, if (u, v) G E then either u ■< v or v ^ ii. We hence may assume that E is a directed tree, 
i.e., E contains only edges (u,v) with u ^ v, in which case we call u the parent of u and u a child 
of u. Note that, except the root, each vertex has a unique parent. We call a vertex a leaf ii it has 
no children. We call G a rooted binary tree if each vertex has at most two children. We distinguish 
binary trees by labelling the children of each vertex as the left child and the right child. 

By 71 we denote the set of finite rooted binary trees with k + 1 leaves, such that the root has 
precisely one child and every other vertex has precisely two children, except for the leaves. Note that 
7o consists of one tree only, which consists of the root, a leaf and an edge going from the root to the 
leaf. Furthermore, put T = UfceN ^ or a ^ ree ™ ^~ we can ^he ver ti ces other than the root and the 
leaves splitting vertices. Note that a tree in 71 has precisely k splitting vertices. For T = (V, E) G T, 
we denote by S the set of its splitting vertices and by L the set of its leaves; hence V = {0} U5UI, 
#5 = k and #L = k + 1. We write T = (0, S, L, E). See Figure [T] for two representatives from T3. 





Figure 1. Two trees in 73. The empty circles represent the leaves L and the full 
circles represent the branching vertices S. 

Now we equip trees with numberings. For k £ No and T = (0, S, L, E) G 71 let /: {0} U S — > 
{0, 1, 2, . . . , k} be a bijection. We call I a monotonous numbering of T if 7(0) = and 7(si) < /(S2) 
for any s±, S2 G S with (si, S2) G E 1 . We extend 7 to L by setting 7(i) = + 1 for any leaf I G L. See 
Figure [2] for an example. The set of monotonous numberings of T is denoted by M{T). 

Now we equip numbered trees with times. For k G No and t > 0, denote by Z^{t) the set of time 
vectors 

t = (t ,...,t k+1 ), where = t < h < ■ ■ ■ < t k < t k+1 = t. (2.1) 

Let us fix a tree T G 71, an ordering I G Af(T) and a time vector t G Z k (t). For 6 = G £7, we 

denote by 

Y^ = {Y r ^:re[t Iiu) ,t Hv) }) 

a continuous-time simple random walk on Z, d with generator re A, starting from zero. We assume that 
the collection (Y^ b,T} )b£E is independent. We consider as the segment of a branching random 
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walk with parent u and child v that arises from a splitting event at time t/( n ), considered until the 
next splitting event at time tjr v y 

Now we compose all these segments of simple random walks according to the tree and define the 
BRW on [0, t] with precisely k splits. Fix the starting site x G 7L d of the branching process. For a leaf 
I G L let = no, . . . , iij-i,Uj = I be the vertices visited by the unique path from to the leaf I, and 
bi = (ui_i,Ui) the corresponding bonds, where j G N. Then we define the continuous-time random 
walk X« = (X?>) reM by 



i-l 



XCO := x+ Y t ^+Y r ^, r G [tj^^tj^i G {!,...,;}. (2.2) 



m=l 



Note that the collection of the random walks (X ( ''); g £ is consistent in the sense that, for any leaves 
I and I', the paths of and X^'^ coincide up to the time tj^) of the vertex u where the tree path 
—¥ I splits from the path — )■ I'; afterwards they are independent given the site X^ = X^ ( _ y The 
separate pieces of the BRW between subsequent splits are denoted by 

X (u, v) = {X^) r&[tl(u)Mv)] = (XV) re[tl{u)Mv)] , (u,v) G E, 

where I G L is any leaf such that the bond (u, v) lies on the unique path from to I. Because of the 
above consistency property of (X^ l) )i e i, the value does not depend on the choice of I. The collection 
of all the path pieces X (u - v) with (u, v) G E is consistent in the sense that X^' v ^ = X^ for any edges 
(it, i?) and (v,u'). See Figure [3] for an example. 

Expectation with respect to the collection (X w )i e i will be denoted by K^' 1 '^. For y G r L d , we 
abbreviate 

*,(T, J, *,»):= f dfE<^[exp( £ £(X<->) dr) ( £ l{X?> = y}], 

(2.3) 

where in the product u is the parent of v. Furthermore, we define 

*, (r, /, t) : = I d?Ef ^ [exp ( f I(v) ^ X r u,v) ) dr ) ( n • ( 2 - 4 ) 
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Figure 3. An example of a BRW corresponding to the monotonously numbered tree 
on the right of Figure [2j 

Finally, we define sequence of numbers Cfc jn for n € N and k = 0, . . . , n — 1 by setting co, n = 1 for 
all n € N and by the recursive relation 

n—k , s 

Ck,n = ^2 ( • ) c k-l,n-U k = 1, . . . ,n - 1. (2.5) 
1=1 

Now we can state the main theorem of this section, which gives us a Feynman-Kac-type formula 
for the functions m n . 

Theorem 2.1. For n £ N and x,y € we /iai>e 

n-l 

m n (t,x) = ^^ c k ^ x {T,I,t), (2.6) 

fc=0 TGTfc l£Af(T) 

and the same formula for m n (t, x, y) with & X (T, I, t) replaced by & X (T, I, t, y). 

Proof. We denote by N(t) the set of particles alive at time t. For a particle u G N(t) let a u and t u 
denote the birth and death time of u, respectively. We put a u (t) = a u At and T u {t) = r u At. If 
u £ iV(t) we write Zg for the position of the unique ancestor of u alive at time s £ [0, t\. If it has no 
children we say that Zi is at the graveyard state, d, for any s > t u . 

Specialising [HR121 Section 2] to our situation, we define a new branching process by imposing 
the following rules: 

(i) We start with one particle at x, which carries n marks (and their positions) 1,2, ... ,n. 

(ii) We think of each of the marks 1,2, ... ,n as a spine and denote by Q the position of the 
whichever particle that carries the mark i at time t. 

(iii) Particles diffuse as under P x , i.e., as independent continuous-time random walks with gener- 
ator At A. 

(iv) A particle at position y carrying j marks branches at rate 2 ] £ > 2{y) and is replaced with two 
new particles. 

(v) At such a branching event of a particle carrying j marks, each mark chooses independently 
and uniformly at random one of the two particles to follow. 
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(vi) Particles not carrying any marks behave as under P x 



We write Qz n) (-) for the corresponding probability measure and Qz n) [-] for the corresponding 
expectation. We call the collection of particles that have carried at least one mark up to time t the 
skeleton at time t and write skel(t). We define D(v) as the number of marks carried by a particle v. 
Let us abbreviate 

rr v (t) , 

(2^)-l)e 2 (ZW)-e (^ ) )) dr' 



A(t) = exp 



E 

ueskel(t) 



a v (t) 



(2.7) 



We now apply the many-to-few lemma [HR121 Lemma 3] for Y = 1 and £ = 1 and obtain 



m n (t,x) 



and 



m n (t,x,y) 



t;6skel(i) 



(2.8) 



Note that the spine trajectory does not undergo a splitting at a branching event, if all the marks 
choose the same child to follow. Hence, we only want to consider splitting events that not all the 
marks choose the same particle to follow. Note that the probability of such event when a particle 
carrying j marks branches is 1 — 2 _ - J+1 . Then the rate of such branching events for particles carrying 



j marks at position y is (2 J — 2)^2 (y)- Accordingly, we define a measure 
(iv) and (v) in the above description of Q x by 



Tin I 



by changing the items 



(iv) A particle at position y carrying j marks branches at rate (2 3 — 2)^2 (y) and is replaced with 
two new particles. 

(v) At a branching event of a particle carrying j marks choose uniformly at random one of the two 
particles to follow conditioned on for each new particle there is at least one mark following it. 

Note that (pH is still valid when Qi"' is replaced by Q^ n) . 

We only prove (|2.6p since the proof of the formula for the moments of m n (t, x, y) is done exactly in 
the same way. We proceed the proof by the method of strong induction. For n = 1, (|2.6p is immediate. 
Now assume that ()2.6[) holds for n replaced by any i € {1, . . . , n— 1}, and we prove that it is also true 
for n. 

We start from the first formula in (j2.8|) and integrate over all values of the time, T, of the first 
branching event under Q x and over all possible branchings. The conditional distribution of T given 
is given by 



{T>t\ C (1) ) = exp ( - |V " 2)£ 2 (C ) ) dr 



(2.9) 



On the event {T > t}, we have skel(t) = {0}, a % (t) = 0, r (<) = t, D(0) = n and Z w = C (1) . Hence, 
we have 

r-t 



A(t)l{T>t} 

(n) 

Integrating this with respect to Q x , we get 



(1) 



exp 



(2.10) 



Q^[A(t)l { T> i} ] 



l.h.s. of (|2~T0|) 



exp / £(A r ) dr 



exp / £(C«) dr 



(2.11) 
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since any spine follows a simple random walk with generator kA. This is the term that corresponds 
to k = in the sum in (12, 6j) . Similarly we can calculate the conditional density of T as 



(T G dti | C (1) ) = exp (- j\i n - 2)6(C^) dr) (2" - 2)6(0 d*i, *i > 0. (2. 



12) 



Let Bi tn _i be the event that at the branching time T, I marks follow the first child of and n — I 
marks follow the second child. Then it is clear that 



n 



1 



I = 1, . . . , n — 1. 



2™ - 2 

So for I = 1, . . . , n — 1, by (12.8)) . we have, using the Markov property at time t\ £ [0, i], 
C ) [^)|^,n-i,T = t 1 ,C (1) ] 

= exp ( T {(2" - 1)6(C W ) - £o(C)} dr) m,(t - h, $>) m n _^ - ti, O 



(2.13) 



(2.14) 



Hence, using (|2.14p . (|2.12p . (|2. 13j) and the independence of the splitting time and the number of 
offsprings, we get 



n-l 



£qJ° A(t)l Bltn _ v T e dtx 



(1) 



n-l 

E 

z=i 



(2.15) 



exp 



£(0 dr) 6(C) m,(t - h, O m n _K* - *i, O dfc 



The induction hypothesis for m/ and m n _; says that 
l-l n-l 

= EE E E E E qh I iCfc ai „- J $ Ct (i)(ri,Ji,t-ti)* ( (i)(r 2 ,i 2 ,t-ti). 

fcl=0fc 2 =0TieT fcl T 2 GTfc 2 hGAfiT!) h£jV(T 2 ) H 

(2.16) 

Let us denote by T (1 ' 2) the tree in Tk 1 +k 2 +i formed by attaching the tree T\ to the left of the unique 
child of the root and the tree T 2 to the right of the unique child of the root. Then the Markov property 
at time t\ gives the following concatenation property of <5: 



E E 

heAf(Ti) heAf(T 2 ) 



exp( / dr)^ w (T 1 ,I 1 ,t-t 1 )^ w (T 2 ,I 2 ,t-t 1 ) 



dti 



Then, integrating both sides of ()2. 15|) with respect to Qj 1 , we get 

q { ; } [A(t)t {Tem} ] = ® { ; } [i-h-s. of cnm] 



n— 1 / \ I— 1 n— i 

n 



z=i 



E , EE E E c ^ c * 



fc 1= fc 2 =0Tier fcl T 2 GT fc; 



2l n-/ E **(r (1 - a) ,i,t). 

/GA/'(T( 1 ' 2 )) 



(2.17) 
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Let T k kl,k kl 1 denote the set of trees in Tk such that the two subtrees of the child of lie in Tk x and 
7ft—ki—l, respectively. By changing the order of the sum we get that the right-hand side of (|2.17p is 
equal to 

n— 1 fc— 1 n-(fc-fei) . . 

EE E E ( jcfci.jcfc-fci-i.n-i E **cw)- ( 2 - 18 ) 

fc=l fc 1= Tg7 -fci,fe-fei-l l=ki+l ^ ' l£jV{T) 

By (12.51) we have 

n-(fc-fci) , . 

E ( 7 ) °ki,l Ck-h-^n-l = Cfe.n- (2.19) 
l=k!+l ^ ' 

This, together with ()2.11j) finishes the proof of (|2.6p . □ 

Remark 2.2. There are also versions of Theorem 12.11 for more general branching mechanisms as 
proposed in Section 11.41 above. Under the additional assumption that X^fceN ^ n £fc converges almost 
surely, one can extend Theorem 12. II to this setting. The main change in (|2.6p is that the terms involving 
the c nj k and £2 must be replaced by a term of the form 

00 

mark V k=2 

where the marks are now taken from a more complex set than {1, . . . , n}. Since the formulas arising are 
much more cumbersome, we abstained from writing them down carefully and proving them. However, 
it is easily seen from the above proof that they have a form which also admits an analysis of the large-t 
limit of the moments of m n in the same way as we do in Section [3j 

3. Proof of the main result 

In this section, we prove the main result of our paper, the moment asymptotics formulated in The- 
orem 11.31 The proof will be crucially based on the Feynman-Kac-type formula for m n given in 
Theorem 1 2 . 1 1 ab ove . Another important ingredient is a large-deviations principle for the local times of 
the branching random walk, which we will provide in Section 13.11 The proof of the lower and upper 
bound of the moment asymptotics are in Sections 13.21 and 13.31 respectively. 

3.1 LDP for the local times of the BRW 

In this section we formulate and prove a large-deviations principle (LDP) for the normalised occupation 
time measures (the local times) of the BRW introduced in Section fl~3| for a fixed tree T = (0, S, L, E) G 
Tk and a fixed monotonous numbering / € Af(T), as the time parameter tends to infinity. 

We define the local times of the BRW as the sum of the times that its random walk segments 
X (u,l,) with (u, v) G E spend in in a given site z G 7L d . More precisely, assume that T G Tk and let a 
time vector t = (to, ... , G Z^{t) be given and define the local time of the BRW in z G 1j d as 

W= E / /W ^)(*)dr. (3.1) 

Then its total mass of is equal to 

™(*):=EW= £(*/(«)-*/(«))• (3-2) 

z&Z d {u,v)eE 
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Hence, we normalise the local times and obtain 

m{t) 

a random element of the set V(Z d ) of all probability measures on Z rf . Fix the starting site x of the 
BRW. Let be the lattice cube of length 2R + 1 centred at x. We consider the periodised local times 

Lf\z):= H z + v)i (3-4) 

ye(2R+l)Z d +x 

a random element of the set ^(T^). Our LDP reads as follows. 

Lemma 3.1. Fix £ No, T € 7/c, I £ M(T) and J?£ff. Furthermore, fix the starting site x of the 
BRW and a vector ? = (so, • • • , Sfe+i) G -^fc(l) and a sequence si — >• s as t — > oo. Then the normalised 
local times Ljg satisfy, as t — > oo, the (full) large deviation principle with scale tm(s) and rate function 
KSft Cl) , where 

SrV)= £ (v^-a/M^)) 2 - (3.5) 

Proof. The special case = is classic and well-known, see [DV75-831 IG77|, IGM98] . Here = (0, 1) 
and m(s) = 1, and the BRW consists of just one random walk with start in x and time interval 
[0, t]. This LDP holds even locally uniformly in s € Zk(l), as is seen from the proof, which uses an 
eigenvalue expansion and the Gartner-Ellis theorem, to turn it into modern notation. This also shows 
that the LDP is the same under the sub-probability measure that conditions on a fixed starting site 
and restricts to a fixed terminal site. 

The general case is an easy consequence of that classical result, as the random walk segments X^ u ' v) 
with (n, v) € E are conditionally independent, after conditioning all the starting site and restricting 
to all the terminating sites, and these are only finitely many. Under this sub-probability measure, 
the normalised local times of each segment X^ u ' v ^ satisfy the LDP with scale t(sj^ — sj^) and the 
rate function in (|3.5h . and Ljg is just an elementary finite convex combination of these independent 
objects. The claimed LDP follows by summing over all the starting and terminating sites, as these 
are only finite sums. □ 



3.2 Proof of the lower bound 

Jensen's inequality gives, for any n,p G N and any x,y £ Z d , 

x)) > (mP{t, x, y)) = (E x (n(t, > (E x (n(t, y)D = (m^(t, x, y)), t> 0. (3.6) 

Now we can apply the result (jl.7p from [GM98] for np instead of p and obtain the lower bound in 
(jl . lip ; see also (|1.12p . Here we recall that the logarithmic asymptotics of (m" p (t, x, y)) are the same 
as for (m™ p (t, x)), as mentioned in subsection 11.21 

3.3 Proof of the upper bound 

Now we give the corresponding upper estimate for the moments (mn(t, x)) (which applies then certainly 
also to (m,n(t, x, y)))- Recall from (|2.6p that (mU(t, x)) is the expectation of the p-th power of the sum 
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of c kin $ x (T,I,t) over k G {0, . . . , n - 1}, T G T k and I G jV(T), where <f> x (TJ,t) is given in (IXlj) . 
Rewriting $ x (T,I,t) using the local times of the BRW, 



z k {t) 



ves 



(3.7) 



it is clear (in the case p = oo, from Assumption ll.2|) that it is only the exponential term involving the 
local times that will turn out to be responsible for the claimed asymptotics, and only the summand 
for k = n — 1 will turn out to give the leading asymptotics. We will dig this term out with repeated 
applications of Jensen's and Holder's inequality. 

We use the inequality 

'J2x t y)<\xrij2(xf), (3-8) 

derived from Jensen's inequality, three times to get 



a — i 



(3.9) 



k=o Te% ieN{T) 

n-l 

= e ° w EE E <*SCW)>- 

k=0Te% leJV{T) 

Hence, the only dependence on t now sits in the last expectation, (<J>§(T, I,t)), and it is enough to 
show that this term satisfies the claimed asymptotics and that it is maximal for k = n — 1. Using 
(|3.8p in integral form for the integral over t and Jensen's inequality for the expectation, we see that 



& x (T,I,t)<\Z k (t)\ 



p-i 



z k (t) 



dt[E£ 



(T,I,t) 



zei d ves 



where = t k /k\ denotes volume of Z k (t). Since this volume term is negligible for the logarithmic 

asymptotics, we only have to concentrate on the integral. Making the change of variables t = ts, we 
see that 



d>£(T,J,i)<e 



o(t) 



ds E 



z k (i) 



ex P (£ e(^))(n6(^) 



(3.10) 



zei 



ves 



At this stage we separate the cases < p < oo and p = oo. 

< p < oo: Our next step is to take expectation with respect to the branching/killing environ- 
ment and to separate the exponential term from the powers of the £2 terms by means of Holder's 
inequality. Fix q, q' > 1 (later chosen in dependence on t) satisfying - + ^ = 1, then we have 



ex P (£ eco^jn 



ze2 



< (E^ [ exp (g £ tizMzj) ] [ n 4 



P/9' 



zG2 



(3.11) 



< (Ef [exp (g £ ^W^))] P ) 1/ '(Ef ■ [ J] )] P ) 



P\ 1/9' 



lies' 
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where we used Holder's inequality twice. Now we show that the second term in the above display is 
negligible. Recall that 151 = k — 1 and that £2 0*0 is i-i-d. in s S Z^. Then using Jensen's inequality 
and Fubini's theorem we get 



E™[n^(^)] p ) w ^ K T ^(n^'(^t)))] i/9 * (3.12) 

where for the last inequality we used the fact that £2 > 0. Now using the inequality x < e x for x > 
we get 

(ef -^'(O)) 1 " < (eP(*-^(0)) W = exp {itf 2 ((fc - 1W)}, (3-13) 

where we recall that H2 denotes the logarithmic moment generating function of £2 (0) . Now we pick 
q = qt = 1 + St and q' t = 1 + 1/ 'st depending on t such that et \ as t — >• 00. In our case where p < 00, 
we have -f^W < p£log£ + O(t) as t — >■ 00, so it is clear that we can choose e t converging to slowly 
enough so that as t — > 00 

jH 2 (l/e t ) 0. (3.14) 

Hence, by H3. 14f) and (|3. 13)) we can conclude that the right-hand side of (I3.12|) is e°W, i.e., the second 
term on the right-hand side of (|3.1ip is negligible. 

Proceeding as in the proof of (jl.7p in [GM98j and using the LDP of Lemma 13. 1) we get that 



exp (q t ^ t{z)tts{z))] P ) llqt < exp (q t ptm{s)) ~2dKptm{s)x(p / n)+o{tj) . (3.15) 



Recall that the LDP in Lemma 13. II holds even uniformly in s£ Z^il) away from zero. Hence, we can 
easily conclude that 

($£(T, J, t)) p < e o(<) / ds exp (-H(q t ptm(s)) - -2dK P tm(s)x(p/ k)) ■ (3.16) 
Jz k (i) y Qt qt f 

By (|1.6|) . for p > we have H(t) S> i as i — > 00 and for p = 0, %(0) = 0. Finally, since qt — > 1 as 
i — >• 00, by Laplace's method we get that, for any T, I and k € {0, . . . , n — 1}, 

<*P(r, /, t)) p < e°W exp (-H(q t pt(k + 1)) - -2dKpi(fc + l) X (p/«)) , (3.17) 

\qt qt > 

as the main contribution comes from s~ = (0, . . . , 0, 1), having total live time m(s) = k + 1. The 
interpretation is that the main contribution to the moments comes from the BRW splitting into k 
particles practically immediately after the beginning. 

Hence, using (|3.9p and the fact that H(t) 3> t for p > and x{p) = for p = once again, by 
Laplace's method ('the largest rate wins') we get 

(m£(t, x)) < e o{t) exp (—H(q t ptn) - —2dKptnx{p/K)) . (3.18) 
\qt qt > 

The proof of the upper bound for < p < 00 is therefore finished by noting that 

exp(-H(gtptnj) = e H(pfn) e° W , t -> 00. (3.19) 
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This is seen as follows. Recall that qt \ 1, and note from [GM98, Remark 1.1(b)] that the convergence 
in (|1.6j) is uniform on [0, 1] and hence also locally uniform on [0, oo). Hence, writing 

1(-H(q t tpn)-H(tpn)) = ~P n( ltH(t) H(pnt) - pnH(t)^ ( q t -lpnH(t) - H(tpn) 



Qt 



(3.20) 



t t 
shows that (|3.19p holds and finishes the proof. 

p = oo: We start from (|3.1U|) . In order to express the p-th power of the expectation, we introduce 



p independent copies , i = 1, . . . , p, of X^'^ and denote by the sum over i G {1, . . . ,p} of 

the local times of these random walks. For z G Z, d define 



and introduce the notation 



G 2 (l,k) 



i=i (u,v)eS 
(e^(0)e 2 (0)*) _ (e«(°)&(0)*> 



^6 (or 



Z,/c G [0, oo), 



(3.21) 



(3.22) 



where the last step used that £2(0) and £o(0) are independent; recall that £ = £2 — £o- From (I3.10P we 
have 

■<fr£(r,/,t)) < e°W / ds(E[exp ( ^ J] &(*) r(z)) 



a o(*) 



a o(t) 



d?E[(exp(^^(^)) IJ6W 



r(z)) 



d ?E[( n 



dsE[ J] (e«(°)^ > W6(0) rW )' 

nG 2 (^(z),r(z)) 



(3.23) 



where we used Fubini's theorem and (|3.22p . 

Note that Zc 1— > G2(l,k) is log-convex for any Z, since it is a moment generating function. Since 
G 2 (/, 0) = 0, it is also easily seen to be log-subadditive. As a consequence, diG?,(l, k) = C?2(Z, k + 1) — 
G 2 (Z, 1)G(Z, A;) is nonnegative, hence Z 1— >■ G2(l,k) is increasing. Hence, since the local times £^{z) 
sum up to tpm(s), we have 

J] G 2 (4| ) ( Z ),r( Z )) < H G 2 {tpm(s),r(z)) < G 2 (tpm(s), £ r(z)) = G 2 (tpro(5), (n - l)p), 

zGZ d zGZ d zGZ d 

since the r(z) sum up to p\S\ = p{n — 1). Note that the right-hand side is < by Assumption 11.21 
Using this fact in (|3.23p , we see that 

(& x (T,I,t)) < e°W / dsE\ TT e H ^W) 
Now, precisely as in the proof of fjl -Tj) in [GM98 , one proves that 



E 



< e 



H (tpm(s))—2dKtpm(s)+o(t) 
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An inspection of the proof, using the uniformity in the LDP in Lemma [3. H shows that this convergence 
is locally uniform in s. Hence, like in the above proof in the case p < oo, we see that Laplace's method 
yields the result, after optimising over s~G Z)~(l) and k £ {0, . . . ,n — 1} (recall that = !)■ 
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